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An integrable spin-ladder model with nearest-neighbor ex- 
changes and biquadratic interactions is proposed. With the 
Bethe ansatz solutions of the model hamiltonian, it is found 
that there are three possible phases in the ground state, i.e., 
a rung-dimerized phase with a spin gap, and two massless 
phases. The possible fixed points of the system and the 
quantum critical behavior at the critical point J = are 
discussed. 

75.10.Jm, 75.30.Kz, 75.40.Cx 



Recently, there has been growing interest in the spin 
ladders for their relevance to some quasi-one-dimensional 
materials, which under hole-doping may show supercon- 
ductivity m . It is well known that the S = 1/2 isotropic 
spin ladders with even number of legs have a spin-liquid 
ground state with an energy gap, while odd-legged lad- 
ders have a gapless spin-liquid ground state. On the other 
hand, generalized ladders including other couplings be- 
yond the nearest-neighbor exchanges, which can inter- 
polate between a variety of systems, can show remark- 
ably rich behavior. |^ Ej In a recent paper 0, Nersesyan 
and Tsvelik predicted a new gapful phase for the two-leg 
spin ladders, i.e., the dimerized phase driven by the four- 
spin interactions, which is essentially different from the 
well known Haldane phase |j) . This observation has been 
demonstrated in a generalized spin-ladder model || by 
constructing the exact ground state. Another interesting 
phenomenon in the ladder systems is the quantum phase 
transition from the gaped phase to the gapless phase, 
which has been studied experimentally in the Heisenberg 
ladder Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 . § 

As is well known, the integrable models provided us 
very good understanding to the correlated many-body 
systems in one dimension. However, a satisfactory in- 
tegrable ladder model, which may play a similar role of 
the Heisenberg chain 0, the one-dimensional (ID) Hub- 
bard model |8| and the supersymmetric t — J model |j] , 
is still absent. The difficulty to construct an integrable 
ladder model is almost the same we encountered in con- 
structing a two-dimensional integrable model due to the 
strict conditions for the integrability. For example, in 
ID, there is only one path connecting two different sites, 
while even for two coupled chains, we have a large number 
of paths connecting two different sites. We note that an 
integrable ladder model with artificial three-spin interac- 
tions has been proposed recently |l(3] and the integrabil- 



ity of a generalized spin ladder without free parameter 
11 was addressed. The latter is more interesting but 
still defys Bethe-ansatz solution. In this letter, we study 
a spin ladder with biquadratic interactions. By properly 
choosing the four-spin coupling constants, we show the 
model is exactly solvable via algebraic Bethe ansatz. The 
model hamiltonian we shall study reads: 
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j=i j=i 

1 N 

+j u i ' ^i+iWi ■ Tj+x) (1) 

3=1 

1 N 

+7 c/ 2^(^ ' TiWi+i ■ T3+1), 

3 = 1 

where aj and Tj are Pauli matrices acting on site j of 
the upper and lower legs, respectively; Ji and J 2 are the 
coupling constants along the legs and the rungs, respec- 
tively; Ux t 2 are the biquadratic coupling constants and 
N denotes the length of the ladder. Without the four- 
spin terms, Eq.(l) represents the ordinary spin- ladder 
model. The new terms in Eq.(l) represent an interchain 
coupling and an interrung coupling, which can be cither 
effectively mediated by spin-phonon interaction or in the 
doped phase generated by the conventional Coulomb re- 
pulsion between the holes moving in the spin correlated 
background as discussed in refs. [2,12]. The importance 
of biquadratic exchange for some properties of Cu0 2 pla- 
quette has been pointed out jl3| and recent experiments 
revealed that such multi-spin-exchange interactions are 
realized in the two-dimensional (2D) solid 3 He fjij ], 2D 
Wigner solid of electrons formed in a Si inversion layer 
plU , and the bec solid 3 He Jl6| |. We note that when 
U 2 = 0, Eq.(l) is reduced to the model considered in 
ref.[2]. For general parameters Ji :2 and U\ t2l the model 
(1) is still non-integrable. However, as we shall show be- 
low, when U\ — J\, U 2 — or U\ = J\, U 2 = —J\/2, 
the model is exactly solvable. We shall study these 
cases through this paper. Not lossing generality, we set 
Ji = U\ = I, J 2 = J and U 2 — U in the following text. 

We study first U — case. This is the simplest in- 
tegrable case but shows the main physics of the system. 
The hamiltonian (1) for U = can be rewritten as 
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In this form, the integrability of the model is still some- 
what hidden. To show it clearly, we note that the 
first term in Eq.(2) can be rewritten as ^2f = i Pj,j+i, 
where Pj,j+i is the permutation operator between two 
nearest rungs. Therefore, the first term of Eq.(2) is 
ST/(4)-mvariant as showed for the spin-orbital model 
Jl7| |. An obvious fact is that Pj,j+i can be expressed 

as P jd+1 = Ea,p X f X i+i> where X f = K >< Pi\ 
are the Hubbard operators and the Dirac states |aj > 
span the Hilbert space of the j'-th rung and are orthog- 
onal (< aj\/3j >= 5 a p). A basic representation of these 
quantum states is |cr|, r? >, where erj, t? = j, J,. However, 
these states are not the eigenstates of the local operator 
<7j ■ Tj . This can be overcome by choosing another basic 
representation 
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U>-U,T>), 
|i>=| T,T>, 



|2>=-^(| t,4>+l U>), 
|3>=|4,4>. 
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The first state denotes a singlet rung and the latter three 
indicate the spin-triplet states of a rung. With these 
notations, Eq.(2) can be rewritten as (up to an irrelevant 
constant) 
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Obviously, the operators N a = Ylj=i A" Q , which denote 
the numbers of the a-rungs in the whole system, are con- 
served quantities. The constant 27 in the last term of 
Eq.(4) indicates a chemical potential applied on N and 
reduces the global SU(A) symmetry of the hamiltonian 
to £7(1) x SU(3). Now we have reduced the hamiltonian 
(1) to an 5[/(4)-invariant spin chain (or equivalently an 
SU(4) t — J model), which can be solved by following 
the standard method There are three branches of 
flavor waves (generalized spin waves) in this system. If 
we choose the reference state as \Sl >= |0i > <&\0 2 > 
<8> • • • <8> |0jv >, these flavor waves describe the spin-triplet 
"excitations". With this reference state, we obtain the 
Bethe-ansatz equations (BAE's): 
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where M x = Ni + N 2 + N 3 , M 2 = N 2 + N 3 and M 3 = 
A3; Xj, Ha and v& represent the rapidities of the flavor 
waves. Note the periodic boundary conditions X^f +1 has 
been used in deriving Eq.(5). The cigencnergy (up to an 
irrelevant constant) of the hamiltonian (4) is given by 
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Obviously, for < 7 < 2, the ground state consists of 
real Xj, Ha and vg closely packed around the origin. That 
means we have three "Fermi seas" and three branches of 
gapless excitations. For 7 > 2, the reference state be- 
comes the true ground state and any flavor excitation 
is gapful. The ground state is a product of the singlet 
rungs, which indicates the dimerization along the rungs. 
The energy gap can be easily deduced from Eq.(6) as 
A = 2(7 — 2). — 2 indicates a quantum critical 
point at which the quantum phase transition from the 
dimerized phase to the gapless phase occurs. At this 
critical point, all the three branches of flavor excita- 
tions are marginal and the low-temperature thermody- 
namics of the system shows non-Fermi-liquid behavior 
as we shall discuss below. For 7 < 0, the singlet rungs 
are unfavorable at low energy scales. For convenience, 
we choose |li > ®|l2 > <8 • • • <8 1 1 jv > as the refer- 
ence state. The BAE's are still given by Eq.(5) but with 
Mi = N a +N 2 + A 3 , M 2 = N 3 +N and M 3 = A . The 
eigenenergy (up to an irrelevant constant) is given by 



=1 X 3 + 4 
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A hidden fact is that there is another critical value 7°. 
For 7 < J1, no singlet rung exists in the ground state 
configuration and the excitations consisting of singlet 
rungs are gapful. In this case, there are only two branches 
of gapless flavor waves and the effective low-energy hamil- 
tonian is equivalent to that of the S'?7(3)-invariant spin 
chain. The ground state consists of two "Fermi seas" (for 
A and fi) with Mi = 2A/3, M 2 = A/3 and M 3 = 0. Wc 
denote the distributions of A and h m the ground state 
as Pi(A) and p 2 (h)i respectively. A singlet excitation can 
be constructed by introducing a v mode and a /Lt-hole Hh 
in the BAE's. We denote further the changes of pi(X) 
and p 2 (h) y ia the v mode and the /z-hole as &p\(X) and 
Sp 2 (p), respectively. From the BAE's (5) we can easily 
obtain 



Spxiu) 



1 



4cosli I - 1 



(8) 



where 8p\(uj) is the Fourier transformation of <5pi(A). 
Combining Eq.(7) and Eq.(8), we derive the minimum 
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energy (corresponding to v — > 0, ///i 
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The critical value Jf. is thus derived from e m i n = as 
J C L = -tt/(4V3) + (Zra3)/4. For J c _ < J < 0, the system 
behaves as for < J < 2. Exactly at the critical point 
J = J1, one branch of the flavor excitations (the sin- 
glet one) is marginal. Therefore, we have three quantum 
phases in these system: A rung-dimerized phase when 
J > J+, a gapless phase with three branches of gapless 
flavor excitations when > J > J1 and another gap- 
less phase with two branches of gapless flavor excitations 
when J < Jf\ We note that the dimerized phase shows 
a long range order 



(10) 



which indicates the condensation of the singlet rungs. 
Under hole-doping, the system behaves as a t — J ladder 
and the singlet rungs serve as Cooper pairs. The mo- 
bility of the singlet rungs under hole-doping may drive 
the system to show superconductivity. Based on the 
above observations, we conclude that J± represent two 
unstable fixed points of the system. In addition, the sta- 
ble fixed points of the system can be conjectured. For 
J > J?, the transverse exchange dominates over the ex- 
change along the legs and the system should flow to a 
fixed point J* = +oo. For J£ < J < Ji, the two unsta- 
ble fixed points J± indicate an intermediate stable fixed 
point Jf. < J* < Jl. For J < Jf. , the singlet excitations 
are eliminated at low energy scales and the system should 
flow to a fixed point J* — — oo, which is equivalent to an 
S'J7(3)-invariant spin chain. The gapless modes in the lat- 
ter phase is mainly due to the high symmetry. Any small 
perturbations of J\ or U\ breaking this symmetry may 
drive it to the Haldane phase as in the S'?7(3)-invariant 
spin-1 chain. 

Based on the BAE's, the thermodynamics of the 
present model can also be derived by following the stan- 
dard method In the gapless phases, the system 
behaves as a Luttinger liquid [pof and nothing is anoma- 
lous. However, at the quantum critical points, the system 
may show non-Fermi-liquid behavior due to the marginal 
excitations. We consider first the zero-temperature mag- 
netic susceptibility for J = J?. Without the external 
field, the ground state is a condensate of singlet rungs. If 
we apply a very weak external field on the system, some 
triplet rungs with S z = 1 appear in the ground-state con- 
figuration while N 2 and N 3 still keep to be zero since the 
levels of these two types of rungs are either lifted (|3 >) 
or unchanged (|2 >). The energy density of the ground 
state in an external magnetic field (H > 0) reads 



E/N= [ A (4-—L T -H)p 1 (\)d\ t 

J-A A ' 4 

where pi(X) satisfies 
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witha„(A) = n/27r[A 2 + (n/2) 2 ] and A 2 = l/(4-H)-l/4. 
For a small H « 1, we have A w and Eq.(12) 

can be solved up to 0{H Z / 2 , A 2 ) as 
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Combining Eq.(13) and Eq.(ll) we readily obtain the 
susceptibility as 
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The low temperature susceptibility and the specific heat 
can also be derived from the so-called thermal Bcthc 
ansatz. fl8|,|l9| Via low-temperature expansion of the 
thermal BAE's we obtain 



C~ Ts, 



(15) 



which indicate a typical quantum critical behavior. 
These results can also be predicted by a simple flavor- 
wave theory with the dispersion relation e(fe) ~ k 2 . We 
note in the gaped phase, the magnetic field can also drive 
a quantum phase transition. At the quantum critical 
point H c = 2(J — 2), similar quantum critical behavior 
can be obtained. 

Now we turn to the U = —1/2 case. The last term 
in Eq.(l) can be rewritten with the basic representation 

Eq.(3) as - T,j[ 2X j° X j+l ~ X j° + 1/4]- U P to an irrel - 
evant constant, we rewrite Eq.(l) as 
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The above hamiltonian deserves similarity to an SU(l\3)- 
supersymmetric t-J model, which still allows Bethe- 
ansatz solution. We choose still |fi > as the reference 
state. The BAE's read: 
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The eigenenergy of the hamiltonian (16) is given by 

Mi 

s = E(^-t + 2 - / - 1 )- 

3=1 X i + 4 



(18) 



The situation is very similar to that of U\ = 1, U = 
case. There are still three phases, i.e., a rung-dimerized 
phase and two gapless phases. For J£ < J < 1/2, 
some triplet rungs are allowed in the ground state. The 
ground-state configuration is described by closely packed 
real v modes and the corresponding A — 3— strings and 
ji — 2— strings: 



Aj = i/ s + i(2-n), n=l,2,3, 



4 



(19) 



For J > J? = 1/2, we get again a dimerized ground state. 
Comparing to the SU(4) case, we find J? is remarkably 
reduced by a negative [/. Notice that a negative £/ in- 
dicates the attraction between two nearest singlet-rungs. 
This attraction enhances the dimerization along the rung 
direction. 

For small positive J or negative J, the triplet rungs 
are more stable than the singlet ones. For convenience, 
we choose |li>®---®|ljv>as the reference state. The 
BAE's lEl read 
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where Mi = N 2 + N 3 + N , M 2 = N 3 + N and M 3 = N Q . 
The eigenenergy takes the same form of Eq.(7) but with 
J — * J — 1/2. J1 can be easily derived from Eq.(20) as 
Jl = 5 — + j In 3. Interestingly, Ji takes a positive 
value in this case. The space of the phase with three 
branches of massless excitations is remarkably depressed 
by the attractive rung-rung interaction. This observation 
strongly indicates that there is a critical point U = U c . 
When U < U c , J± coincide each other and the interme- 
diate fixed point will be eliminated, implying only two 
phases can exist in the system. 

In conclusion, we propose an integrable spin-ladder 
model which exhibits rich physics. This model may play 
a similar role in the spin-ladder systems as the supersym- 
metric t— J model does in the one-dimensional correlated 
electron systems. 
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